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Abstract 

Strong asymptotics of polynomials orthogonal on the unit circle with respect to a 
weight of the form 

m 

W{z)^w{z)Y[\z-ak\^f^'', \z\^l, |afc| = l, /3fc > -1/2, k^l,...,m, 

k=l 

where w{z) > for \z\ — 1 and can be extended as a holomorphic and non-vanishing 
function to an annulus containing the unit circle. The formulas obtained are valid uni- 
formly in the whole complex plane. As a consequence, we obtain some results about the 
distribution of zeros of these polynomials, the behavior of their leading and Verblunsky 
coefhcients, as well as give an alternative proof of the Fisher-Hartwig conjecture about 
the asymptotics of Toeplitz determinants for such type of weights. The main technique 
is the steepest descent analysis of Deift and Zhou, based on the matrix Riemann-Hilbert 
characterization proposed by Fokas, Its and Kitaev. 



1 Introduction and statement of the main results 

Let us first set some notation that will be widely used in what follows. We denote by 
T the unit circle on the complex plane C (circle of radius 1 centered at the origin), and 
O == G C : \z\ < 1} is the open unit disc. If r < 1, let = {z £ C : r < |z| < 1/r}. 
Every oriented Jordan curve or arc 7 induces naturally the notions of left (or and 
right (or "— ") sides of 7. We also denote by a subindex (respectively, "— ") the left 
(respectively, right) boundary values of functions on 7. 

An integrable non-negative function W defined on T is called a weight if 

/ Wiz) \dz\ > 0. (1) 
Jt 
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For each weight W there exists a unique sequence of polynomials (pn (called Szego polynomi- 
als), orthonormal with respect to W, satisfying (pn{z) = Knz"^ + lower degree terms, Kn > 0, 
and 



ipn{z)(pm{z) W{z)\dz\ = 6mn ■ (2) 

We denote by ^n{z) = ipn{z)/Kn the corresponding monic orthogonal polynomials. It is well 
known that they satisfy the Szego recurrence 

<^n+liz) = Z^niz) -C^Kiz), '^oiz) = I , (3) 



where we use the standard notation = z"'<I>„(l/z). The parameters q„ = — <I>n+i(0) 

are called Verblunsky coefficients (also reflection coefficients or Schur parameters) and satisfy 
a„ G D for n = 0, 1, 2, . . . (see QS] for details). 
If the weight W satisfies the condition 

logl^(z)|dz| > -oo, (4) 

then the Szego function of W (see e.g. ^1 Ch. X, §10.2]), 

D{W- z) exp j^^ log W{e^') dO^ , (5) 

can be defined. This function is piecewise analytic and non-vanishing, defined for z ^ T, and 
we will denote by D\ and its values for |2;| < 1 and \z\ > 1, respectively. It is easy to 
verify that 

and in particular, W{z) = \Dc{W; z)\^'^ for z G T. The role of the Szego function in the 
description of the asymptotic behavior of the orthogonal polynomials is well known; for 
instance, 

^.^ ^ ^ D,{W;z) 

n^oo z"^ De{W]00) 

uniformly on each compact set in the exterior of the unit disk. 

In this paper we focus on weights of a specific form. Assume that w is a. strictly positive 
function defined on the unit circle, which can be extended as a holomorphic and non-vanishing 
function to an annulus Ap (0 < p < 1); in fact, a "canonical" form for such a function is 
w{z) = |/(2;)p, z G T, where / is holomorphic and non-vanishing in Ap. For points G T 
and values (3k > —1/2, /c = 1, . . . , m, we define the following weight: 

m 

W{z) = w{z)\{\z-ak?^\ z£T. (7) 

k=l 

It can have zeros or blow up at a^'s, but still conditions and @ are satisfied. We are 
interested in the asymptotic behavior of the corresponding sequences {<I>„(z)}, and {a„} 
when n ^ oo and z G C. 
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The case when ah f3k = (that is, when is a positive analytic weight on T) has been 
studied in where a canonical representation of the corresponding orthogonal polynomials 
in terms of iterates of the Cauchy transform of the scattering function of W, 

S{W;z) = Di{W;z)D,{W;z), (8) 

was derived. Also, detailed asymptotic formulas were obtained. A characterizing feature of 
this case is that the zeros of $n's stay away from T, clustering (with a possible exception of 
a o{n) number of them) at an inner circle determined by the analytic continuation of Dc- 

When any Pk ^ 0,W is no longer a positive and analytic weight on T, and in this situation 
the majority of the zeros of the Szego polynomials are attracted by the unit circle. One of 
the main goals of the paper is to provide asymptotic formulas for $n's valid uniformly on the 
whole complex plane, and in particular, in a neighborhood of a^'s. For partial results in the 
case when all = 1 see |2j. 

In order to state the main results we need to introduce a new piece of notation. 

Let A = {ai, . . . , am}] for the sake of brevity hereafter we call these points generically as 
"singularities" of the weight, although W is regular at a^'s for integer PkS. We fix for what 
follows < 6 < 1 — p, such that additionally S < ^ miuj^j | aj — aj], and denote 

= {z £ C : \z — a^l < S} , = {z £ C : \z — a^l = 6} , k = 1, . . . ,m , (9) 

as well as B = U^^^Bk- Furthermore, given a subset X C C and a value a E C we will use 
the standard notation a ■ X = {ax : x £ X}; consistently, A - X = U'^^^{ak • X). 

Under assumptions on W in ((T)), both Di{W; z) and D(.(W] z) admit an analytic extension 
across T \ we keep the same notation for these analytic continuations. More precisely, 
D\{W; z) is holomorphic in {z G C : \z\ < 1 / p} \ {A ■ [1, and Df,{W; z) is holomorphic 

in {z e C : |2;| > p} \ (^ • (p, 1]); see Section |2l for a detailed discussion. 

First we describe the asymptotic behavior of $„'s away from the singular points of the 
weight: 

Theorem 1 For manic orthogonal polynomials corresponding to the weight W given in 
(O there exist a complete asymptotic expansion 

A:=0 

valid uniformly in C. Each ffe(2;) is a piecewise analytic function, holomorphic in each domain 
specified below. In particular, there exist constants t?^ G T, k = 1, . . . ,m, defined by formula 
below, such that: 

(i) formula 

holds uniformly on every compact subset ofH; 
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(ii) formula 



Di{W;z) ^ Cfc — z , 
/loZffs uniformly on every compact subset of Ap \{B U {A - (p, 
(in) formula 

'''^^"^-^ Z)e(T^;oo) [' + nf^^a, 

holds uniformly on every compact subset o/ C \ B. 

Remark 1 It is well known (see _7^_ as well as Section |21 below) that further terms of the 
expansion (|10|) can be obtained by nested contour integration and the calculus of residues. 
However, the difficulty of the computation increases with k, and we limit our attention to 
the leading nontrivial terms of (|lf 




Corollary 1 For every compact set K C'B there exists N = N{K) G N such that for every 
n> N , every has at most m — 1 zeros on K . 

However, the global behavior of these "spurious" m — 1 zeros can be complicated; we 
describe their limiting set below (Theorem 

In order to formulate the asymptotic behavior of the Szego polynomials in a neighborhood 
of each singular point we define an auxiliary function: if C^/-^ stands for the main branch 
of the square root in C \ {—oo, 0] (that is, C^^"^ > for C > 0), and Ju are the Bessel functions 
of the first kind, then for /3 > —1/2 set 

(-^ 42? re^27r*/3^i/2 (iJ^^^/2(C) + Ji3-i/2{C)) , if C is in the second quadrant, ^^^^ 
\C'^^^ {iJf3+i/2iC) + Jp^i/2{C)) , Otherwise. 

With the notations introduced in Q we have the following result about local behavior of 
the polynomials at the singularities of the weight: 

Theorem 2 Let k G {!,..., m}, \z — ak\ < 6. For z G Bk define 

Cn(.)='-^;iogf-') , (12) 



where we take the principal branch of the logarithm. Then for z £ Bk, 

*"(^' = /l''*'^|^ + , (13) 

where we take the principal branch of the square root in the described neighborhood of z = at- 
The O (1/n) term in H13|) is uniform in the closed disk Bk- 
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n = 24 



n = 51 



n = 75 




n = 101 



n = 128 



n = 150 




Figure 1: Zeros of for several values of n with W{z) = |(z — l){z — a){z — a^)|*^, z G T, 
where a = exp(7ri\/2). Centers of the circles are the zeros of YlT=i Pk'&kO-^'^^ /{o-k — z). 



Remark 2 It will be shown in Section|3]that D^iW] z)TC{Pk] Cn{z)) is a holomorphic function 
in a neighborhood of z = a^. 

Remark 3 In the particular case of /3fc = (a removable singularity) we have 

n{0; C) = C'/' {^Jl{0 + J-i(C)) = y|e< , 

and takes the form 

De{W;oo) V \n 



(cf. Theorem^ (ii)-(iii)). 
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Remark 4 The steepest descent analysis used for the proof of the theorems above ahows to 
find further terms of the asymptotic expansion for the polynomials. However, for the sake of 
simplicity we decided to restrict our attention to the leading terms bearing already non-trivial 
information about the main parameters and the zeros, as it will be shown next. 

For a weight 1^ on T let us define the geometric mean 



1 



27r 



G[W] exp log [W [e^') dO^j . (14) 

If W is given by 0, then G[W] = G[w] = Dl{W] 0) (see formulas below). 

With respect to the Verblunsky and leading coefficients of the orthonormal polynomials 
we have 

Theorem 3 For weight W given in Q there exist constants 1?^ G T, /c = 1, . . . , m, introduced 
in 1)391) below, such that the Verblunsky coefficients a„ defined by the recurrence © satisfy 

«"=-^|;|^+^(;^) ' (15) 

For the leading coefficients Kn of the orthonormal polynomials ipn the following asymptotic 
formula holds: 



7^ fi--E/3i + of^ll > 

G[2Trw] \ ^ V"- / / 



n — > oo. (16) 



Remark 5 Obviously formula ()16|) makes sense for n > X^fcLi/^fci ^^'^ exhibits, at least 
asymptotically, the growing character of the leading coefficients Kj^,. 

A direct consequence of H16() is the asymptotic behavior of the Toeplitz determinants 
related to the weight W. If we define the moments 

dk = (f z-'''W{z)\dz\ , 

then the Toeplitz determinants are 

Vr.{W)'^dethd,.,)l.^ . (17) 



Theorem 4 Under the assumption above there exists a constant x depending on W such 
that 

p„(T^) = X (GpTTu-])" n^"^!'^' (l + o(l)) , n^oo. (18) 
This formula is in accordance with the well known Fisher-Hartwig conjecture (see e.g. |4j). 
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Remark 6 As it was mentioned above, in the case considered, G[w] = G[W], so we may 
replace it; by Vl^ in the right hand side of both formulas (|16() and (|18j) . 

Let us discuss now how Theorems ^^21 reveal the behavior of the zeros of the polynomials 
Qualitatively we can describe the picture as follows: the vast majority of the zeros will 
approach the circle T regularly and radially uniformly along level curves r„ defined below, 
that are asymptotically close to circles centered at the origin with radius 

n-V«. Sin ffularities 

of the weight exert the following influence on the zeros of ^n- those closest to points 
converge to T faster than the rest (their absolute value behaves like c~^^"', with c > 1 
depending on and either leaving a "gap" around (if (3k > 0) or approaching it radially 
(if Pk < 0). Furthermore, a bounded number of "spurious" zeros may wander inside the unit 
disk; each will have at most m — 1 of these zeros (Corollary and their global behavior 
will depend in particular on the relative positions of points a^. 
Let 



n 



^ (jk^kClk 



k = l " 



(19) 



where S{W;z) is the scattering function defined in These curves are well defined: al- 
though S{W;z) is multivalued in a neighborhood of T, by (pTTj) below, |5(VF;z)| is single- 
valued and positive in Ap. Furthermore, given an analytic function / and e > denote 



Z{f) {z : fiz) = 0} , Zeif) '^{zGC: nun^ \t-z\<e} , (20) 



and let 



Tnie) = i z £ C : min It — zl < e 
' ter„ 



be the e-neighborhood of r„. 
Let also 



?i = 1, and 6*^ = — (argafc - argai) , k = 2,...,m, (21) 
zvr 



and 



m 



^^^^^ def y Mk_ ^2.^{n+l)e, (22) 

Theorem 5 There exists < e < 5 such that for all sufficiently large n G N every "pie slice" 
of the form 

z e r„(e) : ^ < arg(z) < ^ | C \ (B U Z,(7^„)) , 

(n) (n) 

with appropriately chosen a £ M, contains exactly k = k2 — ki zeros of z| , . . . , , 
satisfying 

|.f^| = l-i^ + ofiV (23) 
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arg(.ffi)-arg(.f)) = ^ + o(-l). (24) 
See Figure |21 for an illustration of the statements of Theorem El 



n = 15 n = 30 n = 45 




-1 1-1 1-1 1 



n = 60 n = 75 n = 90 




-1 1-1 1-1 1 



Figure 2: Zeros of for n = 15,30, ... ,90 with W{z) = \{z - l){z - a){z - a?)\'^, z G T, 
where a = exp(7ri\/2). Points 1, a and are indicated with asterisks. For comparison with 
the prediction of Theorem [S] we plot in each case the level curve r„ defined in H19() . 

We can also be more precise about the accumulation set of zeros of <5n's, 

fc>ln>fc 

on compact sets of B. Assertion (i) oi Theorem ^ shows that the structure of ZnB depends 
on the relative positions of ai, . . . , am on T. Without loss of generality we may assume that 
^1 = 1, ^2; • • • ! is the maximal subset of . . . , 9m} linearly independent over the rational 
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numbers Q. Then there exist unique r^j € Q, k = 1, . . . ,m, j = 1 



■ . . ,v, such that 



V 



0k = '^rkjOj, k = l,...,m 



(obviously, r^j = 5kj, for 1 < k < v). Let qi = 1, and m > 2, < pk < qk £ Z,, k = 2, . . . ,m 
be such that 



Theorem 6 Let teZnB. If v = 1 (that is, if all 6^ G QJ, then there exist < Sk < qk, 
Sfc G Z, k = 1, . . . ,m, such that 



Moreover, every point in ID that belongs to the manifold given by solutions of l\27^ when 
X2, . . . vary in M, is an accumulation point of the zeros o/{<^„}. 

Corollary 2 (i) If v = \, then Z n B is a discrete set of a finite number of points. 

(a) If V = 2, then Z 01} is an algebraic curve of degree < m. In particular, if v = m = 2, 
then Z n D is either a circular arc (if \f3i\ 7^ \l32\), or a diameter in D formed by the 
perpendicular bisector of the segment joining ai with 02 (if \f3i\ = \P2\)- 

(Hi) If V > 2, then Z H D is a two-dimensional domain bounded by algebraic curves. 

Remark 7 From the method of proof it follows in fact that the zeros of the orthogonal 
polynomials will be equidistributed on the manifolds described above. 

Compare these statements with the numerical results depicted in Figure 01 Take note that 
we plot all the zeros of for n = 1, . . . , 150, in order to reveal the structure of Z inside D. 

Finally, in order to describe the behavior of the zeros that are closest to the singularity 
of the weight we must consider again function TC{P;C) introduced in (fTTl) : let us denote by 
h{(3) its zero of smallest absolute value such that Re(/i(/3)) > and Im(/i(/3)) > 0. According 
to (IHSJ, z = —h{f3) is also a zero of 7i{(3; (). 

Theorem 7 Among the zeros of ^n{z), the closest to the singularity ak of the weight are 
those, given asymptotically by 



Tki = — mod Z, k = 2, . . . ,m . 

Qk 




(26) 



If V ^"^j then additionally there exist X2, . . . , X^, G M such that 




(27) 



z = Ofc e' 



,2ih(h)/n 



(1 + 0(1)) and z = a^e 



2ih(Pk)/n 



(1 + 0(1)) 



approaching symmetrically with respect to the radius • (0, 1). 
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1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 
-1 

-1 -0.5 0.5 1 

Figure 3: Zeros of for n = 1, 2, . . . , 150 with W{z) = \z — l\^/^\z — a|~^/^, z G T, where 
a = exp(7rz\/2). Points 1 and a are indicated with asterisks. Observe different features of the 
behavior of the zeros described in the text: the zero of the weight at z = 1 repels the zeros of 
the polynomial, while the singularity at z = a attracts a zero of ^niz) that approaches z = a 
radially (cf . Remark |^ . Also the arc of the circle along which the zeros of $n remaining 
inside B cluster is clearly visible (cf. Corollary Ej) • 

Remark 8 Using formulas (9.6.3) of J] we see that C = —ih{(3) must satisfy the following 
equation: 

Vi/2(C)- Vi/2(C) = 0, Re(C)>0, Im(C)>0, (28) 

where Ip is the modified Bessel function. Some facts about the zeros of the function in the 
left hand side of (|28|) have been kindly provided to us by M. Muldoon ^1]. For instance, for 
/? > this function is strictly negative in (0, +oo) (see (1^) and thus has no real positive 
zeros. Furthermore, for — 1/2 < /? < this function has apparently a unique real positive 
zero, which is a monotonically decreasing function of [3. Complementing these results with 
numerical experiments (see e.g. Remark^] in Section we claim that: 

• for /? > 0, Re(/i(/?)) > and \h{[5)\ > vr; 

• for -1/2 < (3 <0, Re(/i(/3)) = and \h{P)\ < vr. 

In particular, for —1/2 < (3 < 0, z^ = z~, lying (at least, asymptotically) on the radius 
afe-(0,l). 
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-0.3 -0.25 -0.2 0.85 0.9 0.95 1 1.05 1.1 



Figure 4: Zeros of $150 (dots) in a neighborhood of z = a (left) and z = 1 (right) for 
VF(2;) = \z — ll^^'^lz — a|~^/^, z G T, where a = exp(7ri\/2). Points 1 and a are indicated 
with asterisks. Centers of the circles are the zeros of 'H(— 1/3; Ci5o(-z)) and 7^(1/6; Ci5o(-z)), 
respectively. 

Remark 9 We can compare the rate of convergence of the zeros of to the unit circle. 
While, according to Theorem 121 the bulk of the zeros tend to T with a speed whose leading 
term is n"^/" (Theorem ISJ , the zeros closest to the singularity are attracted by T with a 
rate, roughly speaking, equal to c"^/"", where c = e^^^'^^^''^ > 1. Taking into account also 
that \z^ — Ofcl = ^\h{(3)\ (1 + o(l)). Remark |H1 and (HU, it explains the influence that a^'s 
exert on the neighboring zeros of both pushing them to T and "repelling" them from (if 
/3 > 0) or attracting them to (if /? < 0) the singularity (cf. Figure 01 and zoom in Figure j^. 

Remark 10 The method of proof of the statements above can actually handle a more general 
situation, when the weight of orthogonality depends on n, 

m 

Wn{z)= w{z) Yl\z - akl'^f^"'" , zeT, (ofc G T for = 1,... ,m), 
k=l 

assuming that for values f3k{n) > —1/2, k = 1, . . . ,m, the limits 

lim = /3a: > , k = 1, . . . ,m , 

n—>oo n 

exist. 

The structure of the rest of the paper is as follows. In the next section we study in 
more detail the properties of the Szego and the scattering functions of W. The Riemann- 
Hilbert characterization and the nonlinear steepest descent methods of Deift and Zhou (see 
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e.g. IBIEI) as well as the monograph |HI), are carried out in Sectional which allows to prove 
the announced results in Section 0J 



2 Szego and scattering function for W 

Following the notation introduced above, let Di{w]z) and Dc{w;z) be, respectively, the in- 
terior and exterior values of the Szego function defined by © for the analytic component w 
in 0. Since w is holomorphic and non-vanishing in Ap, then both Di(w; z) and De{w, z) ad- 
mit a holomorphic extension across T, and maintaining the same notation for these analytic 
continuations we have 

Di{w;z) 

—— r = W(Z) . 

This formula gives in a certain way a canonical analytic extension of w from T to the annulus 
Ap. 

It is convenient to construct explicitly the Szego function for the modified weight W] with 
this purpose we select the singe-valued branches of the corresponding functions as follows: 
the generalized polynomial 



q{z)'^J[{z-a,)M^ 



(29) 



fe=i 



is a single-valued analytic function in C \ (U^j^Ofc • [l,+oo)), for which we fix the value of 
Q'(O). Consequently, q{l/z) is single-valued and analytic in C \ (U^-^a^ • [0, 1]), and 

1 



q{l/-z) 



9(0) 



If we consider each ray ■ [1, +00) oriented towards infinity, then 

g+(z) = e-^^^'=g_(z), z G • [1, +00) , k = l,...,m. (30) 
In the same fashion, if ■ [0, 1] also have the natural orientation from the origin to a^, then 



q{l/z) 



q{l/z) , zGak- (0,1), k = l,...,m. 



With this convention we can write the Szego functions for the modified weight W: 

^2^-^ D,{w;z) 



Di{W;z) 



q [z) 

qm 



Di{w;z), D,{W;z) 



(0) q{l/-z) 



2 • 



(31) 



(32) 



Then Di{W; z) is holomorphic in Di/^ \ (U^^^a/c • [1, l/p)), Dc{W] z) is holomorphic in {z G 
C:\z\ >p}\(U^^iafe-(/9,l]), and 



A(W^;0) D,{w-Q) 



Dc{W;oo) = Dc{w;oo) > 0. 



(33) 
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Using the definition in Q we see that r is related with the geometric mean defined in (|14|1 
by 

r-2 = G[W] = G[w] . (34) 

Furthermore, formula 

Di{W;z) 



W{z) = q\z)q{l/-z) w{z) = -^Jy^^ (35) 

is valid and provides an analytic continuation of the weight W to the cut annulus Ap \ 
(U^j^Ofc • {p,l/p)). By (|Hn|) - (|3H) . with the orientation of the cuts toward infinity we have 



[A(t^;^)]+ = e-2'^*'^MA(VF;z)]„ , z e • (1, 1/p) , fc = l,...,m, 
[De{W;z)]^ = e~^^'^>'[D,iW;z)]_ , zGafc-(p,l), k=l,...,m. 
Recall that the scattering function for w, 

S{w; z) = Di{w; z)Dc{w; z) , 
is holomorphic in the annulus Ap. By 



(36) 



5 ( w; - ) = — , for z G Ap , \S(w, z)! = 1 on T . 

V zj S{w;z) 

With the definition Q and formulas (|32)) we have 



SiW;z)=DiiW;z)D,{W;z) 



q{z) 




^q^{0)q{l/-z)^ 

that is also analytic and single- valued in the cut annulus Ap \ (U^-^a^ • (p, l/p)). By 

S+{W;z)=e-^^'f'^S.{W;z), z e ■ {p,l/ p) \ {ak} , k = l,...,m. (37) 

It is straightforward also to check that S{W; z) is bounded in a neighborhood of each a^. 

Recall the definition in Q . With our assumptions on 5 > the closed disks Bk U Ck are 
disjoint, and S{w; z) is analytic in each B^- Furthermore, motivated by (|57|) we define 

c s dcf /e''*^'=5(H^;z) , if z G and arg(z) > arg(afe) , 

SkiWiz) = < .„ k=l,...,m. (38) 

\e-^*^*5(VF;z), if z G and arg(z) < arg(afc) , 

Then Sk{W;z) is holomorphic in Bk, k = 1, . . . ,m, and we define 

^k=SkiW;ak)£T, k=l,...,m. (39) 
These constants have been used in the formulation of several results in the previous section. 

Remark 11 Formally 

" n f 5(u;;afc), 



but the selection of the right branch in each case should be specified. So the definition by 
(|39|) is preferred in order to avoid ambiguity. 
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3 Riemann-Hilbert analysis for orthogonal polynomials 



We assume the unit circle T oriented counterclockwise. The starting point of all the analysis 
is the fact that under assumptions above conditions Q can be rewritten in terms of a non- 
hermitian orthogonality for and : 

^n{z)z'^-^-'^ dz = 0, for A; = 0, 1, . . . , n - 1 , 



W{z) 



dz 



0, A; = 0, l,...,n-2, 

k = n — 1 . 



By standard arguments (see e.g. [3] or |S], as well as the seminal paper ^ where the Riemann- 
Hilbert approach to orthogonal polynomials started), 



Yiz) 



r ^n{t)W{t)dt \ 
Im Jj t^{t - z) 

Kn-l f ipl.,{t)W{t)dt 



(40) 



\ ^' ' I h t^{t-z) J 

is a unique solution of the following Riemann-Hilbert problem: Y is holomorphic in C \ T, 



y+(z)=y_(z) 



1 W{z)/z' 
1 



zeT. 



where I is the 2x2 identity matrix, and 

O 

Yiz 



O 



1 \z-ak\^'^^ 
1 \z-ak\^P^^ 
1 l\ 



1 1 



lim Y{z) 



if < ^ 
if A >0; 









as z — > Ofc, z G C \ T, and A: = 1, . . . , m. 

In order to perform the steepest descent analysis as described in |H] (see also ^U]) we build 
a series of explicit and reversible steps in order to arrive at an equivalent problem, which is 

solvable, at least in an asymptotic sense. We will use the following notation: 173 







-1 



is the Pauli matrix, and for any non-zero x and integer m, x 



0-3 



X 
1/x 



3.1 Global analysis 

If we define 



Hiz) 



z^"'''^, if \z\ > 1, 
/, if \z\ < 1, 



(41) 
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and put T{z) = Y{z) H{z), then T becomes holomorphic in C \ T, with limz^ooT{z) = /, 
and satisfying the jump condition 

on T. Obviously, T exhibits the same local behavior at A as Y. 

Let 7i be a closed Jordan contour, piecewise analytic, entirely contained in the cut annulus 
{z £ C : p < \z\ < 1} \ U^j^Ofc • (/>, 1), except for points a^: 7; n T = ^. Let also 
7e = {l/z : z G 7i} the contour symmetric to 7i with respect to T. We take both 71 
and 7e oriented counterclockwise. Let be the connected component of C \ 7i containing 
the origin, and fioo the corresponding unbounded component of C \ 7e. Furthermore, we 
denote by 

Oi = D\n^, = {l/z : z G Jli} 

(see Fig. El. 




Figure 5: Opening lenses. 
Then W is correctly defined in both ^\ and fie, and we can set 



Kiz^ 



1 

^z'^/Wiz) 1 J 

1 0^ 
l/{z'^W{z)) 1 



if z G ilo U floo, 
, if 2 G Oi, 

, if z G ilo- 



(42) 
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Using K we make a new transformation: U{z) = T{z)K{z). Matrix valued function U is 
holomorphic in C \ (T U 7i U 7c), lim^^oo U (z) = I, and 

U+iz) = U.iz)Juiz), zGTU7iU7e, 



where 

W{z 
-1/W{z) 

1 
^z^'/Wiz) ij ' 

1 
l/iz'^Wiz)) 1 

Moreover, the local behavior for U (z) as z ^ at from fig U Vt\ is now 



Ju{z) 



if z e T \ ^, 
if z G 7i \ ^, 
if z G 7c \ ^. 



(43) 



U{z) 



O 

o 



1 k-afc|2AA 
1 |^_a^|2A.J 

\z-ak\~'^l^^ I' 
|z-afc|-2A 1 



if A <0, 
if A > . 



With De{W; •) and Di{W; ■) defined in (jSH) and constant r introduced in (jHSJ, let 
r/I)c(TV;z)^'^= 



N{z) 





-t/A(VF;z) 



A(VF;z)/t\ /A(W^;2 







r 

-1 0, 



if > 1 , 



if Ul < 1. 



(44) 



By has the same jump on T as U , and by (1321), it exhibits the same behavior at 

infinity. Hence, U{z)N~^{z) tends to / as z — > oo, is holomorphic in C\ (TU7eU7i), and has 
jumps across these curves asymptotically close to /, except when we approach the singular 
set A. We have to handle the behavior at each individual singular point by means of the 
local analysis. 

3.2 Local analysis 

Let us pick a singular point G A. For the sake of brevity along this subsection we use the 
following shortcuts for the notation: a = a^, 

/3 1^' /?^, B = Bi, C = Ce (where Be, Ce and 6 were 
defined in Q), and B^ = {z £ B : arg(2:) > arg(a)}, B^ = {z £ B : arg(z) < arg(a)}. We 
also write = J7j H B, where j G {i, e, 0, oo}, and analogous notation for curves: 
etc. Finally, S{W; z) = St{W] z) (see the definition in (IHHll). 

Our goal is to build a matrix P{z) == P{a, (5; z) meeting the following requirements: 

(PI) P is holomorphic in B \ (T U 7i U 7c) and satisfies across T U 7i U % the jump relation 
P+{t) = P^{t)Ju{t), with Ju defined in (gnj. 
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Figure 6: Local analysis in B. 

As a first step we reduce the problem to the one with constant jumps. Let us denote 
Li = a • (1 — 5, 1) and Le == a • (1, 1 + 5), oriented both from z = a to infinity (see Fig. 
Let w^^'^{z) and z^/^ denote the principal holomorphic branches of these functions in b, and 
(cf. (ESI)), 

W^/\z) = q{z)'^{l/^w^/\z) 

with q defined in ^I^ . Then W^^"^ is holomorphic in B \ a ■ {1 — 5, 1 + 6) , and according to 
©-(Ell), 

}^-.-^'^ „„r,. and '^=e-'^ onF.. 
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Thus, if we define 



\{f3-z) 



-Trip 



-7ri/3^1/2(- 



z]z 



n/2 



ZjZ 



-n/2 



ZjZ 



-n/2 



z e (OeUOoo)n^+, 
z G (OiUOo)ne+, 



(45) 



and set 



Riz) = P{z) X{f3; zy^ , z G 6 \ (Fi U Te U T U 7i U 7e) , 



(46) 



we get for R the fohowing problem: R is holomorphic in 6\ (riUrcUTU7iU7e), and satisfies 
the jump relation R-\-{z) = R-{z)Jr{z), with 



if z G T, 



Jr{z) 




if zG (7ine+)u(7ene-), 
if zG (7ine-)u(7ene+), 

if z G Ti U Te . 

Moreover, R has the following local behavior as 2 ^ o: if /? > 0, 



Riz) 



O 



o 




if z G U il, 



if z G u r^i. 



and if /? < 0, then 



R{z) = O 



As it could be expected, this problem resembles very much the one we face during the 
local analysis for the generalized Jacobi weight on the real line (see flHl Theorem 4.2], and 
jl2j). We take advantage of the results proved therein in order to abbreviate the exposition. 

Let us define Sj {xe*'^/^ : x > 0}, i = 1, . . . , 8; we take Sj oriented towards the origin 
for i = 3,4,5, and the rest towards infinity (see Fig. EJ. This splits the plane into eight 
regions, as indicated there, marked with Roman numbers / to VIII. Using Hankel functions 

(i) 

H{, and modified Bessel functions Ky and ly we build a piece-wise analytic matrix- valued 
function = ^'(/?; •),/?> —1/2, in the following manner (where C}^"^ denotes the main branch 
in (—00, 0]): 
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S3 

\ /// 

\ 

\ 

IV ^\ 



II 



7( 



V 



/ 

/ 

/ VI 

S5 



N VIII 

\ 

\ 

\ 

\ 



VII 



\/ 



E7 



For C e I, 



Figure 7: Local parametrix. 



/3^ 



e-{/3+3)™3^ 



For C G II, 



*(/3;C) 



e 2 



e 2 



e 2 



For C e III, 



M/(/3;C) 



62 



For C G IV, 



r(i) 



(2) 



e(/3+3)™3_ 
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For C e V, 



/-<.(-c) -<4(-c)^ 



(1) 



-(/3+i)7ri<T3 



(51) 



For C G VI, 

m, c) = 

For C G VII, 

*(/3;C) = 

And finally, for C G VIII, 



e 2 



(3'KicFZ 



(52) 



(53) 



(2) 



,(/3+i)7ri<73 



(54) 



Proposition 1 ([18j, Theorem 4.2) Function ^' = ^'(/?; •) defined above is holomorphic 
in C\ IJ^=i ^i; ^^^^ exhibits the following jumps: 



^+(C) = ^-(C) 



1 

-1 



, /orCGS4US8, 



1 

e"2-^/3 1 



/or C G Si U 



*+(C) = ^-(C) (^2^/3 J), /orCGSaUSv 



*+(C) = ^-(C) 
Additionally, if (3 > then as C — > 0, 



e 



'Trip 



for C G S2 U Se. 



(55) 
(56) 
(57) 
(58) 



O 

o 



and if P < then 



z/ C is in domains II, III, VI, VII, 
otherwise, 
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Consider in C \ (—00, 0) the transformation 

f{z) = -i log{z/a) , 

where we take the principal branch of the logarithm. Then / is a conformal 1-1 map of B 
onto a neighborhood of the origin. Moreover, T is mapped onto M oriented positively, and 
we may use the freedom in the selection of the contours deforming them in such a way that 
/(7i) and /(7c) follow the rays Sj in B with odd indices i (dashed lines, see Figs. Inland [71). 
By construction, matrix 

matches the jumps and the local behavior of R in B. Since a left multiplication by a holo- 
morphic function has no influence on the jumps, and taking into account ()46() . we see that 
matrix P can be built of the form 

P(z) = E{z) ^ ^ fiz)) A(/3; z)-^ , (59) 

where E is any holomorphic function in B. An adequate selection of E is motivated by the 
matching requirement P{z)N~^{z) = I + 0(n~^) on the boundary C. 
For the sake of brevity let us denote 

C = ^f{z) = -i^log{z/a), (60) 

(we omit the explicit reference to the dependence of C from z, a and n in the notation). 
Matching condition can be rewritten as 

N{z)X{P;zr^{(3;Cr\ zeC. (61) 

The key idea is to replace ^'(/3;C) by its leading asymptotic term as C — > 00. The complete 
expansion at infinity of the entries of is well known (see e.g. Chapter 9]), so we 

can insert it in formulas ()47 |) - H54() . This computation has been carried out in jl8l Section 
4.3]; we can formulate the result therein by defining the matrix- valued function Q: 

g(f3- e'?^3e"^'^''3e"3^™^ if C is in the first quadrant, (62) 
= eX'^3g-iC^3g|/37ri<73^ ^£ ^ ^Yie second quadrant, (63) 

if ( is in the third quadrant, (64) 

if C, is in the fourth quadrant. (65) 

It is easy to check that G is holomorphic in each quadrant and matches the jumps of ^ on 
E2, S4, Eg and Eg (solid lines in Figurel?!), given in and 



E{z) 



I + Oi - 

n 
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Lemma 1 ([18j) For function 4'(/3;C) defined in (|i7|) - (f5l|) . 



*(/?;C) 



1 / 1 



V2 V-^ 1 



OO .1, 



a(/3;C), (66) 



as C ~^ CO, uniformly for in eac/i quadrant. Here, the constants sp^k o-nd t^^k (^Te given by 



Sf},k 



(67) 



where 



(4zy2 - l)(4jy2 _ 9) . . . (4^,2 _ - 1) 



/n particular, 



V2 



-i 1 



if) 
2C 



/3 i 
z -15 J 



+ 



g(/3;C), C^oo, 



(68) 



uniformly for Q in each sector. 

Taking into account (|6T|) and (|68|) . it is reasonable to set in (|59|) 



i?(z)^^'iV(z)A(/3;z) 



1 / 1 



\/2 \-i 1 



a (AC) 



(69) 



Proposition 2 Matrix valued function E defined in 1)69^ is holomorphic in B and has there 
the following representation: 



E{z) 



S{W-z) 



■ la 



J_ ( i 1 

V2 V-1 



(70) 



where we take the principal branch of the square root. 



Proof. The verification of formula (|70j) is straightforward (and analyticity of S is a direct 
consequence of it). Under transformation (|6U() we have the following correspondence: 

z £ (Qi U Jlo) n 5^ 44> C is in the first quadrant, 

z £ {Qi U r^o) n 44> C is in the second quadrant, 

z £ (Jlc U r^oo) n 0^ 44> C is in the third quadrant, 

z G ($7e U r^oo) n 5^ 44> C is in the fourth quadrant. 

Assume for instance z € (ili U l^o) H in which case, according to (jlH) and (|45j) . 

Ar(^) = f , A(/3; z) = e--'^W'/\z)z~^/^ , 
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and by 

Gathering these elements in (|69|1 we arrive at (|7r)|l . The analysis in the rest of the quadrants 
is similar. 

□ 



Corollary 3 Matrix P{z) = P{a,f3; z), 



P{a,(3;z) 



def 



I f i 1 



V2 V-1 



-0-3 



(71) 



with C given by solves the Riemann-Hilbert problem (P1)-(P2) defined at the beginning 

of this subsection. 

Thus, it remains to check the matching condition (P3). 

Proposition 3 LetP{z) = P{a,f3;z) be given by (fTT]). Then for z G C, function P{z)N^^{z) 
has the following asymptotic expansion: 



P{z)N-\z) =1 
00 .h. 



Sf3,k 



{-lfr~^a^S{W-z)tp^k' 



(72) 



^^2^^K' \r-a--[s{W-z)y\,,, (- 
where constants sp^k cind tp^k were defined in (jH7|) . In particular, 

( 



P(.)Ar-i(z)=/ + ^ 



15 -r-2a"5(H^;z)' 

1 

-/3 



so that P{z)N-^{z) = l + 0{n-^) for z eC. 
Proof By JSH), 

Piz)N~Hz) =E{z)^{l3-C) [N{z)mzr] 
By definition of E in (|69() we have 

N{z)\{(3;zy^ = E{z 

Thus, 



(73) 



3 1-1 



P{z)N-\z)=E{z)^{l3-0 



i?(z)-i. 



(74) 
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By m, 



J_ / 1 -i 

V2 [-i 1 



^(/3;C) 
[I + V{0]g{P;C) 



V2 [-^ 1 
_]_ ( I -i 

J_ ( I -i 

V2 V-^ 1 



a (AC) 

a (AC) 



where V is given by the asymptotic expansion 

v(n V ( ^-'^^''^''^ -'*^A 

k=l ^ VH-lj r/3,fc Sf^^k J 

and constants s^^^ and t^^^ were defined in ()67p. Using it in ()74|1 and taking into account 
(|7n|) we get 





f 1 - 




V-i 1 


T3/2 




' 1^ 


( 


.-1 0, 



-1 



and formula ((7^ follows. Furthermore, since = 2/?^, and t/j^i = 2/3, we obtain ^TS\i . 



□ 



3.3 Final transformation 

With the notation introduced in Q and with P{a, (3; z) defined by (|71j) let us take 
P{z)=P{ak,Pk]z) for zG^fc\(TU7eU7i), A; = l,...,m, 

and put 

. det fc/(z)iV-n^), for z e C \ (i? U T U 7e U 7i), 
\?7(z)P-i(z), for z E 5\ (TU7cU7i). 

It is easy to show that this transformation is well defined, since the inverses exist. Matrix 
S is holomorphic in the whole plane cut along 7 U C, where 

7=(7cU7i)\i? and C = A 
(see Fig. |H|, S{z) I as z ^ 00, and if we orient all C^'s clockwise, S+{t) = S-(t)Js, with 



(P{z)N-^{z), 

1 0^ 
tV(z'^5(VF;z)) 1^ 
1 -z"5(VF;z)/t2^ 
1 



if z G C7, 
if z e 7e \ 

if z G 7i \ 
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Figure 8: Jumps of S. 



It is clear that the ofF-diagonal terms of Js on 7; \ i? and 7e \ i? decay exponentially fast. 
On the other hand, by (f75)) . Js{z) = I + 0{l/n) for z G C. So the conclusion is that the 
jump matrix Js = I + 0{l/n) uniformly for z G 7 U C. In fact, (|7^ gives us the complete 
asymptotic expansion of Jg in negative powers of n. Then arguments such as in [HI El El 
(see e.g. Section 7.2]) allow to show that S itself has an asymptotic expansion in negative 
powers of n. The main observation is that if we define 



where we integrate along contours 7 U C with the orientation shown in Fig. El then M defines 
a bounded linear operator acting in ^^(7 U C) (with respect to the Lebesgue measure), with 
the operator norm 



where "const" depends on 7 U C only (see e.g. (HI Chapter 7] for details). In particular, for 
all sufficiently large n, operator (1 — M) is invertible, and we have the following 

Proposition 4 Let n he such that (1 — M) is invertible, and denote 





(76) 



M||£,2_^^2 < const|| J5 — I\\]^2 




(77) 
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Since for all n sufficiently large, ||M||j;^2^j;^2 < 1, the auxiliary function can be computed 
in terms of the Neumann series, 

oo 

SO that 

oo 

S = Y, S^^^ , with 5(0) = / , 5(^+1) = K (^i^')) , i e N . (78) 

j=0 

In other words, S can be recovered from J5 by nested contour integration. 

Second observation is that, taking into account the exponential decay of the off-diagonal 
terms of Js onjc\B and 7i\-B, we can restrict the integration in (|76j) to contours C, replacing 
then Js by PN~^; in this way, only exponentially small terms are neglected. Pluggin the 
asymptotic expansion of PA^~^, given by (f72|) . into equation (|78|) allows to find successively 
the terms s^. In particular, for terms S^^^ in H78|) we have 

sB(.)=o(i-) 

locally uniformly in C. Furthermore, 

00 /X 



k=l 



uniformly in z, where s^s are piece- wise analytic functions in C\C Let us determine explicitly 
the first nontrivial term Si. By (|76() and H78|) . 

^nw N x^/.N 1 / (P(t)N-^t) - I)dt 

S'^^'{z) = K (/) = - — : * -^-^^^ — h exponentially smah terms, 

where integrals are taken counterclockwise. Hence, we have 
Corollary 4 Matrix S satisfies 

k=i •^'-fe ^ 
locally uniformly for z G C \ (7 U C), where integrals along Ck 's are taken counterclockwise. 
Define foi z G B^D Ck {k = 1, . . . ,m) 

/ -r-^alSkiW;z) 

^k{z) = : 1—. r 

nlog(z/a,) \^^2„-n5~l(^.^) 

SO that by 

P{z)N-Hz) - I = Tk{z) + O (J^y zeCk. 
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Since J^k{z) is a meromorphic function in BkUCk with a simple pole at z = Ofc, by the residue 
theorem, for z ^ B^U Ck, 

2711 Jq^ t — Z t=ak t — Z 

where 



Fk{z) = 

Consequently, 

1 / {P{t)N-^{t)-I)dt fFfc(z) + iizeCXWk, 



2vri Jc, t-z [ {Tk + + ^(l/n^), if z G ^fe . 

In particular, in (|79j) . 

Now we are ready for the asymptotic analysis of the original matrix Y (and in particular, 
of its entries (1, 1) and (2, 1)), that we perform in the next section. 

4 Asymptotic analysis 

Unraveling our transformations we have 

y,.\S{z)N{z)K-\z)H-\z), ifz€C\B, 

\s{z)P{z)K'\z)H-^{z), ifzeB. ^ ^ 

Proof of Theorem Q and Corollary Q Using (|79j) in (|81j) we readily obtain (|10j) . In order to 
prove the rest of the statements we must analyze the consequences of (jHJ) in each domain 
(see Fig. El . 

In \ -B we have (cf. ^ and (jSI)) 

Hence, Y{z) = S{z)N{z), so that 

^nW = -^;(^5i2(-)> ^2i(.) = -^^5.2(.). 
Taking into account (|8()|) and recalling that = Yn we obtain 

.„(.)._L^i(f«i.o(i)). ..fi„\B, 

rA(VF;z) n ak - z \nj J 
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Figure 9: Domains for the asymptotic analysis. 



which yields formula (i) in Theorem ^ 

Observe that the leading term in the right hand side is a rational function with at most 
m — 1 zeros. In consequence, for all sufficiently large n each ^niz) can have at most ni — 1 
zeros on a compact subset of the unit disk (Corollary^. 

Let us analyze the asymptotic behavior of the polynomials close to the inner boundary 
of the unit circle, but still away from the singular points A. For z G $7i \ i? we have 




Hence, 



Y 



(z) = S{z) 



z''Di{W;z)/{TW{z)) Di{W;z)/T 
-T/Di{W;z) 




z^D,{W;z)/t Di{W;z)/T 
-T/Di{W;z) 



) 



Analyzing Yn we obtain that 




(83) 



which coincides with the formula of statement (ii) of Theorem 
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The exterior asymptotic is analyzed likewise. For z G fig \ -B we have (cf. (|IT|) . (|l2|) and 

m) 

N{z) = {{D,{W;z)/tY\ = (-1/(z"V(z)) ?)' ^"'(^)=^"'''- 

Unraveling the transformation we conclude that formula ()83() is valid locally uniformly also 
in J7e \ B, that is, holds in a neighborhood of T away from the singular set A. 
For z G $7oo \ B we have 

iV(z) = ((i?e(H^;z)/r)'^% K = /, i/-i(z) = z"'^3. 

Hence, ^(z) = S{z) {{z'^Dc,{W; z)/t)^^ , and statement (iii) of Theorem^is a direct conse- 
quence of formula (|5n|) . □ 

Proof of Theorem\^ Now we turn to the asymptotic analysis for <I>„'s in a neighborhood of 
a singularity a^, where we have to use the expression 

Y{z) = S{z)P{z)K-\z)H-\z) 

with P{z) = P{ak,f3k; z) for z £ Bk- 

Assume that z G fli D B/. and arg(z) > arg(afc) in such a way that C, = —i^ log (z/ak) 
belongs to sector I in Figured Then ^ is given by (|47)) . which along with (|45|) replaced in 
(|7Hl yields 

<T3/2 

0-3 



\ (T3/Z 

^^^^^ m^J M(C) (e-/^W^-V2(,),"/: 



where 



M(C) 



(C) 



Thus, 



1/2 



1/2 



^'12(^) 



2^/2 



^2 



1/2 



''^ ^ 2V2' Ufc(H^;z) ' ^" i ^ ' 
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Taking into account that K and H are as in (|82|l we see that 

{P{z)K-\z))n = Pii{z) + Pi2{z) 



W{z) 

le'^ ( '-^al^\'' (Mn(C) +^M,.(C)) , 



{P{z)K'\z))2l = P2l{z) + P22{Z) 



z" 



W{z) 

o \ 1/2 

2 'Ti \ ' 



ic'^' ^ L . «fc"T7^ (*M2i(C) +M22(C)) 



Thus, 



yii(z) = Sii{z){P{z)K~\z))i^ + 5i2(z)(P(z)i^-H^))2i 
= (P(z)i^-i(z))n + 
By [H formulas 9.1.3-9.1.4], 

Mn(C) + ^Mi2(C) = 2 (iJ^^+i + J^^„i) (C) , 
so that for z in the domain considered 

where C, = —if log(z/aA;), and 7i{(3;C) has been defined in (fTTI) . Function 7i{(3]C) is holo- 
morphic in the complex plane cut along the positive imaginary axis iM+, and by ^ formula 
9.1.35], 

W+ C) = e^-^ W_ (/?; C) , C e zM+ , (85) 

if the positive imaginary axis is oriented towards the origin. 
Taking into account that 

/ SkiW;z) \ ^ ie^'^'^/^D.iW; z) , ii z e Bk and arg(z) > arg(afc) , ^ ^ ^ 

\^ W{z) j |e-^*/3fe/2De(VF; z), ifzeBk and arg(z) < arg(afc) , ^■■■^rn, 

we see that (fTT?|) is valid for z £ QiCi Bk and arg(z) > arg(afc). 

Since transformation z i— > maps the ray ■ (0, 1) onto the positive imaginary axis 
oriented towards the origin, from formulas (|36|) and (|85|) it follows that D(.(W; z)7i{Pk', Cn) is 
single- valued in a neighborhood oi z = ak- Hence, by uniqueness of the analytic continuation 
we must conclude that (fT3|) is valid in fact in the whole Bk- □ 
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Proof of Theorem\^ For the Verblunsky coefficients we have also 



Or, 



1 / 1 \ 

fc=l ^ ^ 



n — > oo 



and it remains to take into account that £ IR, £ T, and G T in order to arrive at 
formula (fT5|) . 

Furthermore, observe from (|4(H) that the leading coefficients k„ of the orthonormal poly- 
nomials (fn have the following representation in terms of matrix Y: 



1 



Since in Oq \ -B we have Y{z) = S(z)N{z), with 

N{z) = 

then 



Y2iiz) 



A(W^;^)/r 
r/A(W^;^) 

S22{z) , 



so that 



_2 



n-l - ^ '522(0) 



and formula (|16j) follows from (|34)) and (|80)) . 



□ 



Proof of Theorem^ One of the main connections of Toeplitz determinants T>niW) defined 
in (|17|) with the orthogonal polynomials on the unit circle is that they can be expressed in 
terms of the leading coefficients Kn using the following formula (see e.g. |15| Theorem 1.5.11]): 



Vn-i{w) kI ■ 

Since Vq{W) = do = <f z~'^W{z)\dz\, we obtain that 



-2 
'3 



Fix iV e N such that N > J2T=i formula (fTH|) . for n> N. 

Vn{w) = do n '^f n 



do n 

j=N 



T 

2^ 



n -1 
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where we denote £i = do Y[f=i^ '^j ^- Hence, 

^ ^ i=N I ^ k=l ^ 



£ 



27r 
72 



71 n 



n 

j=N 



j + 



1^^ 

k=l 



-1 



n 

j=N 



1 + 



o in 



This last infinite product is convergent to a constant, which we denote by £2- Hence, 

-1 



VniW)=£i£2 



2tt 

72 



n -j^E/'i 

j=N L •' k=l 



2vr\" nn + 2)T{N + l-YJ=iPl 



r(iv + i)r(n + 2-Er=i/3D' 

Gathering all the constants in x and using Stirling formula for the asymptotics of the Gamma 
function we obtain that 

_j n^^=,Pl (1 + 0(1)) , n^oo. 



It remains to use formula (|34() in order to set the proof of Theorem 0] 



□ 



Proof of Theorem\^ Let e > 0. We may rewrite formula (in) of Theorem ^ using the 
notation (j^ : 



aJ+'D,(M'-;0) 



z"5(W^;z)-i7^„(z)[>+F„(z). 



(87) 



If z G r„(e), then there exists a constant C = C{e) such that 

l^n(^)|<-^. 

We conclude the proof using standard arguments involving Rouche's theorem (see e.g. proof 
of Theorem 4 in [T^). □ 

Proof of Theorem\^ Let K he a compact set. For n large enough and z G if, we have 
|z|" < n-2, so that by 

znKcf]\J z{iZn), 

k>ln>k 

with TZn defined in ((211). Thus, it is sufficient to describe all the possible limit points of {Gn}- 
Observe that with the notation introduced in Section ^ before the formulation of Theorem 

ini 



IZniz) = > exp (27ri(n + l)Ok) = 7 exp 27r« > r^j (n — m + 1) 



k=\ 



Z- Ok 



k=l 



Z-Qk 
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Let V = 1; this means that all 9k G Q, and 9^ = r^i = Pk/Qk niod Z, k = 2, . . . ,m. Then 
clearly all possible limits of TZn{z) are given by equation l^^ . 

Assume v > 2. By Kronecker's theorem (also known as Kronecker-Weyl theorem, see 
e.g. Ch. Ill]), since 9i = 1,92, ■ ■ ■ ,9j] are rationally independent then (and in fact, if and 
only if) sequence 



jieN 



is dense and uniformly distributed in the {v — l)-dimensional torus T'' . In particular, for 
any real numbers X2, ■ . . , Xy there exists a sub-sequence A C N such that 

lim e^'^*^^" = e^'^*^^ , j = 2,...,v. (88) 

neA 

In fact, we can say more: 

Lemma 2 Let 9i = 1,92, ■■■ ,9^ be rationally independent, and let rkj € Q, j = 2,. . . ,v, 
k = 1, . . . ,m. Then for any real numbers X2, ■ ■ ■ , X^ there exists a sub-sequence A C N such 
that 

lim e^^^^*^^^" = e^^^'-'^^^^ j = 2,...,v, k = l,...,m. (89) 
Indeed, for rkj £ Z this statement follows trivially from (|88|) . If l/rkj £ Z, then 



and we can specify the single-valued branch of the r^^^-th root in the neighborhood of 
exp(2TTiXj) in such a way that the corresponding limit in ((HHl) holds. Combining these 
two observations we obtain Lemma |21 which shows that the set of limit points of {TZn} is 
given by the left hand side of (|27|) . Now the statement follows for v > 2. □ 

Proof of Corollary\^ There is a finite number of numbers Sk S [0, Qk) H Z; for each possible 
combination of s^'s, the left hand side in (|26|) is a rational function with denominator of 
degree <m — 1. Now the statement (i) follows. 

If f = 2, then Z n D is a manifold parameterized by a continuous parameter X2 £ M, 
which shows that it is a curve. It is easy to check that its equation is a polynomial in two 
real variables of degree < m. Furthermore, \i v = m = 2, then equation ()27l) is equivalent to 



f3i'di{z-a2) ^ ^2^iX2 
P2'&2{ai - z) 



I3i'di{z - 02) 



f32'&2{ai - z) 



which reduces to |z — ai| = |/?2| |z — 02!- This proves (ii). 

Finally, for v > 2, the set Z H D is a manifold parameterized by at least two continuous 
parameter Xi S M, showing that generically it is a two-dimensional domain. Its boundary is 
again an algebraic curve or a union of algebraic curves. □ 

Proof of Theorem^ It is a straightforward consequence of the asymptotic formula ^1'6\ . 
that shows that in the neighborhood Bk of Ok the zeros of match (up to a 0{l/n) term) 
those of Ti.{Pk] Cn)- D 
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Remark 12 We don't know any explicit formula for h{P), although it can be easily computed 
numerically. In order to find a good initial value for we can use the continued fraction 
expansion [li Formula 9.1.73]: 

Mz) ^ 1 

•^-i(^) 2u/z ^' 

2{u + l)/z 

In particular, truncating at the second term and equating to zero we can take 

^ V(2/? + 3)(6/3 + l) + i(2/3 + 3) 



as a reasonably good approximation for any iterative (say, Newton-type) zero-finding method 
of computation of h{P). As an illustration, in the table below we compare the values of /io(/3) 
and h{p) for /3 = -1/4 and /3 = 1, . . . , 5: 



p 


-0.25 


1 


2 


3 


4 


5 




0.69f 


2.96 + 2.5i 


4.77 + 3.5? 


6.54 + 4.5i 


8.3 + 5.5« 


10.04 + 6.5i 


h{(i) 


0.68i 


3.73 + IMi 


5.08 + 0.87i 


6.34 + 0.79i 


7.56 + 0.74i 


8.75 + 0.71i 
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